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Abstract 

New type of nonsingular oscillating solutions for the Universe described by cos- 
mological equations of gauge theories of gravity is presented. Advantages of these 
solutions with respect to existing nonsingular solutions within framework of general 
relativity and gauge gravity are discussed. It is shown in particular that these solu- 
tions have nonzero measure and stable on contraction stage unlike usual nonsingular 
solutions. 
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1 Introduction 



In conventional cosmological scenario the Universe has a beginning in time, so 
called Big Bang, associated with initial singularity of Friedmann equations [1]. 
Singular state with infinite values of energy density and curvature is unaccept- 
able from physical point of view and possibly marks the edge of applicability 
of General Relativity (GR). In spite of some progress in last decades, attempts 
to resolve this problem within GR had no success [2], [3]. 

Correct description of the early Universe when the energy density was very 
high should be given by particle physics [2]. Unfortunately experimental and 
observational data are still insufficient to discriminate between different mod- 
els of elementary particles physics. However in most models scalar fields play 
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an important role. One of the main discovery in modern cosmology is a pos- 
sibility to explain a number of puzzles of Big Bang model by inclusion of 
short stage of exponential expansion of the Universe, so called inflation [2], [4]. 
During this stage the energy density of the Universe is dominated by scalar 
field (or fields) and transition to usual hot Universe is achieved by creation of 
relativistic particles at the end of infiation [5]. 

Unfortunately even in inflationary models there is little chance to avoid cos- 
mological singularity. The evolution of the Universe is usually considered from 
the planckian moment when description in terms of classical theory of gravita- 
tion becomes possible. At contraction stage the singular state is unavoidable 
in most models [3], [6]. Some nonsingular sohitions exist [6]- [10] but they are 
not of great interest since their measure is equal to zero and they cannot be 
used to construct viable nonsingular cosmological models. 

Alternative way to resolve problems of Big Bang model could be so called os- 
cillating model [11] where the Universe never reaches singularity. In the recent 
literature such models are debated in connection in particular to ekpyrotic 
[12] and cyclic Universe models [13]. 

From the other hand, deviations from Riemannian geometry can occur at 
very high energy density. In the case of more general spacetime GR should be 
replaced by gauge theories of gravity (GTG) that are natural generalizations 
of GR [14] . It was shown that nonsingular cosmological solutions within GTG 
exist in the case of spinning matter, usual matter [15] and matter in the form 
of scalar field [16]. Moreover, when energy density drops enough the deviation 
from GR becomes negligible and at late stages the Universe can be described 
by Priedmann equations that possess to generalize standard cosmology and 
resolve the problem of initial singularity. 

In this paper we present new type of nonsingular cosmological solutions for 
the Universe filled by nonlinear scalar field within gauge theories of gravity. 
We show that our solutions have some advantages with respect to previously 
known nonsingular solutions. The paper is organized as follows. In section 2 
cosmological equations of GTG are given. In section 3 oscillating nonsingular 
solutions are presented. In section 4 stability of our solutions is discussed. In 
section 5 discussion is given. The last section contains conclusions. 



2 Cosmological equations 

Homogeneous isotropic models in GTG are described by generalized Fried- 
mann cosmological equations [17] that were obtained first within framework 
of Poincare gauge theory of gravity. Same equations are valid in metric-affine 
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gravity. They read 
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where R{t) is the scale factor in the Robertson- Walker metric, k = —1,0, +1 
for open, flat and closed models respectively, Mp is planckian mass, p{t) and 
p{t) are energy density and pressure, (3 is indefinite parameter, and a dot 
denotes the differentiation with respect to time ^ . 

These equations contain the only indefinite parameter (3 that is a combination 
of gravitational lagrangian coefficients, having a dimension of inverse energy 
density. Solutions of cosmological equations have essentially non-einsteinian 

behavior when the energy density exceeds the critical value I3~^. Experimental 
constraints on parameter /3 are very poor and one can say only that /3~^ < M^. 

Besides equations (1-2) gravitational equations of GTG lead to the following 
relation for torsion function S and nonmetricity function Q 

s-lQ--~W-Pip-m- (3) 



In Poincare gauge theory of gravity Q — and equation (3) determines the 
torsion function. In metric-affine theory of gravity equation (3) may describe 
three kinds of models: in the Riemann-Cartan spacetime Q = 0, in the Weyl 
spacetime 5* = 0, in the Weyl-Cartan spacetime {S ^ 0,Q ^ 0), the function 
S is proportional to the function Q [18]. 

The conservation law has the usual form 

p + 3H{p + p) = 0, ^(^) = |' (4) 

where H{t) is a Hubble parameter. 

In the case of nonlinear scalar field minimally coupled to gravity the energy 
density and pressure are 

p^l^'it) + Vi^), p^l^\t)-V{^), (5) 

^ Hereinafter the system of units with ^ = c = 1 is used. 
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where V{(fi) is effective potential of the scalar field ip{t). 



Then taking into account (4-5) the cosmological equations (1-2) can be rewrit- 
ten in the following way [9], [16], [19] 
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where Z{t) = 1 + /? [p\t) - W{^{t))i V = ™, V" = 

The conservation law (4) in the case under consideration is reduced to the 
scalar field equation 

(p + 3Hip = -V. (8) 



3 Oscillating solutions 



Numerical investigation of nonsingular solutions of equations (6-8) was carried 
out in [9], [16], [19], [20]. It was shown in particular that it is possible to obtain 
nonsingular solutions if parameter (3 is negative and potential energy of the 
scalar field dominates its kinetic energy at the bounce. In the following we shall 
assume simple effective potential usually considered in inflationary models [2] 

V{^) = \\^\ (9) 



where A the a coupling constant. 

The constant A must satisfy the condition A ~ 10~^^ [2] to obtain density 
fluctuations after inflation as small 

as ^ ... 10-4. 
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Figure 1. Oscillating solution with negative parameter /3, corresponding to closed 
cosmological model. 

The typical solution in the case of closed model is shown at fig.l. Initial con- 
ditions for this solution are A = 10~^^, tpo = ■ 10~^ Mp, ipQ = 0, (3 = 
—10^^ Mp'^. The duration of expansion is 5 ■ 10^ 

Initially the scalar field is placed at some distance from the value = and 
the kinetic part of the energy density is zero. Therefore at the beginning the 
effective equation of state is p ~ —p that is one of requirements for inflation. 
At the same time the Hubble parameter is zero and its derivative is positive 
that means transition from contraction to expansion. In course of evolution 
the scalar field goes to zero; potential energy is transformed into kinetic one. 
The Universe expands. 

When the scalar field reaches the minimum of effective potential V{lp) its 
energy is equal to the kinetic energy. The expansion is replaced by contraction 
in that moment. Then the scalar field slows down and the energy density is 
transferred back to the effective potential and we arrive to the initial state. The 
function (3) oscillates with the same frequency and sign as Hubble parameter. 

Oscillations of the scalar field are usual for inflationary models in the frame- 
work of GR, but the stage of oscillations has nothing to do with inflation itself. 
It corresponds to the stage of reheating after inflation during which usual par- 
ticles are created from the scalar field [2]. In our case however oscillations in 
scalar field occur together with oscillations of the whole Universe. 

Note, that in the model under consideration, condition of inflation H <^ is 
satisfied during expansion. The duration of the phase of expansion, however, 
is close to the inverse Hubble parameter, i.e. to the age of the Universe. Thus, 
the expansion rate is very small to allow one to speak about the real inflation. 

It is clear, that usual mechanism of reheating cannot work in this model be- 
cause it was developed for monotonically varying background. However, re- 
cently a very efficient mechanism of preheating was found [21], when the only 
one oscillation is necessary for particles to be created. If some mechanism 
produce scalar field at the stage of contraction, the solution will describe the 
regular cosmological model in which very rapid transition from contraction to 
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expansion appears when the energy density of the Universe is concentrated in 
the scalar field. The rest of the time the Universe is usual Friedmann one. 

It is necessary to note that particular form of the effective potential (9) is not 
so important. As numerical calculations show oscillating solutions are quite 
generic feature for GTG. The only restrictions are that there must be a local 
minimum of the effective potential and it should be sufficiently fiat. 



4 Stability with respect to initial conditions 

In this section we show using intuitive approach that our solutions have one 
advantage with respect to known solutions obtained in GR as well as GTG, 
namely our solutions are stable on the stage of contraction. 

With this goal we vary initial conditions, namely values of (/^o? 0o and Hq at the 
beginning (starting at some moment on contraction stage) and see what hap- 
pens with solutions. Example is given at fig. 2 for finite time intervals^ . Both 
solutions are bounded, so singularity never appears. They are most sensitive 
to variation of initial value of the derivative tpo but if condition ^ ^(v'o) 
is satisfied solutions look very similar. Sensitivity to initial value of the Hub- 
ble parameter is very weak. Note that the frequency of oscillations depends 
on initial conditions even for small variations of (fo. 




Figure 2. Cosmological parameters in oscillating solutions with different initial con- 
ditions. Both solutions are bounded that is an indication of stability of oscillating 
solutions. 

In general case together with our oscillating solutions singular solutions are 
present. However ranges for initial conditions allowing bounded solutions are 
finite although quite narrow (at least in ipQ direction). These finite regions 
represent 'islands of stability' were nonsingular solutions exist. 

Finally we plot 3-dimensional phase trajectories for our solutions with varying 

^ It can be seen that the amplitude of oscillations of cosmological parameters is 
decreasing in time. If one follows solutions for a longer time it can be found that 
the amlitude itself oscillates reminding modulations. The frequency of modulations 
is typically 10^ times smaller than the amplitude of oscillations. 
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initial conditions. It is clear that phase trajectories lay at some 2- dimensional 
hypersurface (probably slice). Each trajectory is a closed curve that means 
that solutions are indeed oscillating. These closed curves do not remind circles 
and therefore oscillations are not harmonic and likely cannot be described by 
analytical functions. 




Figure 3. 3-dimensional phase portrait of oscilating solutions. They all lay at some 
2-dimensional curved hypersurface (or slice). 

5 Discussion 

In this section we would like to carry out comparison to existing nonsingular 
solutions obtained previously within framework of GR as well as GTG for 
cosmological models with single scalar field. 
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Nonsingular cosmological solutions for the Universe filled by scalar field within 
GR are known since a long time [7] and are discussed in recent paper [9], [22]. 
It was shown that such solutions exist in GR only for closed models. Unfortu- 
nately their measure is equal to zero and moreover on the stage of contraction 
they are unstable [6], [23]. For most solutions in course of contraction effec- 
tive equation of state becomes p = p and this leads to singularity. Therefore 
such solutions are of certain interest in connection with the initial singularity 
problem but they do not cure for final singularity. 

Moreover, on the stage of expansion there is transition from vacuum dominated 
expansion regime to Friedmann expansion due to particle creation process. It 
is well known that such process is irreversible [23]. Appearing and subsequent 
domination of the scalar field on contraction stage is an obscure question. 

Nonsingular infiationary solutions with nonlinear scalar field were obtained 
within framework of GTG first in [16] in the case of initial conditions used 
in chaotic infiation [24]. In subsequent papers [9], [20] such solutions were dis- 
cussed to some extent. The crucial property of such solutions is presence of 
inflationary stage during expansion; at the same time quasi-de- Sitter contrac- 
tion stage is present in all of them. 

It is necessary to stress that similar solutions exist within the framework of 
GR [10] most of them are singular on contraction stage since equation of state 
p — p is much more probable. Careful investigation of this problem for the 
case of flat cosmological models with massive scalar held in GTG [19] lead 
to similar conclusions ^ . One can perform simple numerical analysis similar 
to ours and change values of (^q, (po or Hq on contraction stage. It can be 
shown that even httle changes lead to unbounded solutions. This is a typical 
property of nonlinear equations. Therefore in our opinion oscillating solutions 
have important advantage although there is no inflation on expansion stage. 

Finally it should be stated that our solutions can be obtained in the wide 
regions of initial conditions. Moreover, since maximal value of energy density 
depends on parameter /3, oscillating solutions can be obtained for the case 
when p <^ Mp, the Universe can never enter epoch of quantum gravity and 
can be described by classical theory of gravitation. 



6 Conclusions 

In this paper we give a new type of cosmological solutions of GTG for the 
Universe filled by nonlinear scalar field. Numerical analysis shows that these 

^ However effective equation of state is p = | in this case. 
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solutions are stable on contraction stage unlike solutions with inflationary 
stage. If any mechanism of reheating can work in our case, it can be possible 
to build nonsingular cosmological model on the basis of our solutions. 
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